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On two-dimensional percolation clusters at the percolation threshold, we study {a{MB,r)), the 
average conductance of the backbone, defined by two points separated by Euclidean distance r, of 
mass Mb- We find that with increasing Mb and for fixed r, {a{MB,r)) asymptotically decreases to 
a constant, in contrast with the behavior of homogeneous sytems and non-random fractals (such as 
the Sierpinski gasket) in which conductance increases with increasing A'Ib- We explain this behavior 
by studying the distribution of shortest paths between the two points on clusters with a given Mb ■ 
We also study the dependence of conductance on Mb slightly above the percolation threshold. 

PACS numbers: 5.45.Df, 64.60.Ak, 64.60.Fr 
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I. INTRODUCTION 



There has been considerable study of the bond perco- 
lation cluster considered as a random-resistor network, 
with each occupied bond having unit resistance and non- 
occupied bonds having infinite resistance [pl^p| ■ In two di- 
mensions, the configuration studied is typically an L x L 
lattice and the conductance is measured between two op- 
posite sides which are assumed to have infinite conduc- 
tance M-M ■ The backbone of the cluster is then defined 
as the set of bonds that are connected to the two sides 
having infinite conductance through paths that have no 
common bond. 

At the percolation threshold, the backbone mass scales 
as (Mb) ~ L'^'' with de = 1.6432 ± 0.0008 and in 
this "bus bar" geometry is strongly correlated with L. 
The total conductance of the backbone as a function of 
L has been studied extensively and has been found to 
scale as {a) - i"'^ with /t = 0.9826 ± 0.0008 0. 

Recently, the distribution of masses of backbones de- 
fined by two points, i.e., backbones defined as the set 
of those bonds that are connected by paths having no 
common bonds to two points separated by distance r 
within an L X L lattice, has been studied m^ . One finds 
that when r <^ L, there is a very broad distribution of 
backbone masses for a given r. Figure ffl illustrates some 
typical percolation clusters and their backbones defined 
in this configuration. Because of the broad distribution 
of backbone masses we have the opportunity to study 
the conductance between these two points separated by 
a fixed distance, r, as a function of the mass of the back- 



bone defined by these points. 

One might expect that, for fixed r, the average con- 
ductance would increase with increasing backbone mass 
because there could be more paths through which current 
can fiow. In fact, we find that the average conductance 
decreases monotonically with increasing backbone size, in 
contrast with the behavior of homogeneous systems and 
non-random fractals in which conductance increases with 
increasing Mb- We explain our finding by noting that the 
conductance is strongly correlated with the shortest path 
between the two points, and then studying the distribu- 
tion of shortest paths between the two points for a given 
Mb- This analysis extends recent studies of the distri- 
bution of shortest paths where no restriction on Mb is 
placed &|2|. 



II. SIMULATIONS 

Our system is a two-dimensional square lattice of side 
L = 1000 with points A and B defined as ^ = (L — 
r/2, 500), B = [L + r/2, 500). For each realization 
of bond percolation on this lattice, if there is a path of 
connected bonds between A and B, we calculate (i) the 
length of the shortest path between A and B, (ii) the 
size of the backbone defined by A and B and (iii) the to- 
tal conductance between A and B- We obtain data from 
100,000realizationsfor eachof 8 valuesof r (1, 2, 4, 8, 16, 
32, 64, and 128) at the percolation threshold, Pc = 0.5. 
We bin these results based on the value of the backbone 
mass, Mb by combining results for all realizations with 
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2" < Mb < 2"+^ and taking the center of each bin as 
the value of Mb ■ 

In Fig. H(a), we plot the simulation results for the av- 
erage conductance {a{MB,r)) and find that the conduc- 
tance, in fact, decreases with increasing AIb- The de- 
crease is seen more clearly in Fig. 0(b), in which we plot 
scaled values as discussed below. 



backbone. Figure || illustrates how the size of the back- 
bone constrains the possible values of the shortest path. 
For all values of Mb , a section of each plot in Fig. |j(a) 
exhibits power law behavior. In Fig. ^(b), we show the 
distributions P(£|AfB,r) for different r and a given Mb- 
In Fig. Q(c) we see that when scaled with r''min the plots 
collapse, so we can write P{i\MB, r) in the scaling form 



III. SIERPINSKI GASKET 



P{£\MB,r) 



(1) 



In non-fractal systems, the conductance increases as 
the mass of the conductor increases. We next consider 
the average conductance on the Sierpinski gasket, a non- 
random fractal, the first three generations of which are 
illustrated in Fig. ||(a)-(c). Because the Sierpinksi gas- 
ket is not translationally invariant, the analogue of the 
average conductivity between two points in the percola- 
tion cluster is the conductivity averaged over all pairs 
of points separated by distance r . At each successive 
generation, there are two types of pairs: (i) pairs which 
correspond to pairs in the previous generation (e.g., A 
and B) and (ii) pairs which do not correspond to pairs in 
the previous generation (e.g., D and E). It is obvious that 
as we move from one generation to the next, the conduc- 
tance between pairs of type (i) increases because there 
are more paths between the points then in the previous 
generation. On the other hand, the conductance between 
the pairs of type (ii) are lower on average than between 
the pairs present in the previous generation because on 
average the shortest path between the two points is longer 
than between the pairs in the previous generation. How- 
ever, for any given r, the shortest path between any two 
points has a fixed upper bound independent of the gen- 
eration. Due to this bound on the shortest path, the net 
result is that the average conductivity increases with suc- 
ceeding generations. This is shown in Fig. 0(d) in which 
we plot the average conductivity calculated exactly for 
generations 1 to 6 for r = 1, 2 and 4. 



An expression for ^ can be found by recognizing that we 
can write the well-studied distribution P(i\r), the prob- 
ability that the shortest path between two points sepa- 
rated by Euclidean distance r is £, independent of Mb, 
as 



/•oo 

P(i\r)= / P{e\MB,r)P{MB\r)dME 



(2) 



where (i) P{MB\r) is the distribution of backbone masses 
given distance r between the points which determine the 
backbone and (ii) c^ is the lower cutoff on Mb given L 
P{MB\r) has the form |l| 



m.k)-if^ 



[r « L] (3) 



where cLb is the fractal dimension of the backbone and 
Tg = d/ (Ib is the exponent for the blob size distribution 
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P(£|r) 



'91 



(4) 



where dmin is the fractal dimension of the shortest path. 



Since t ~ r'^min and Mb ' 
the lower cutoff q in Eq. 



^^ , implying £ 
scales as 



M^™'" 



eds/d^ 



/dE 



(5) 



IV. SHORTEST PATH DISTRIBUTION 



In order to understand why the average conductance of 
the percolation backbone decreases with increasing Mb, 
we must (i) recognize that the conductance is strongly 
correlated with the shortest path between the two points 
and (ii) study P{£\MB,r), the distribution of shortest 
paths between the two points for a given backbone mass. 
Hence we next create the P{£\MB,r) probability distri- 
bution, binning the data logarithmically by taking the 
average over samples centered at log2 i- 

Figure 0(a) shows the results of the simulations for 
P{i\MB,r) for r = 1 for various backbone masses. The 
plots collapse, the only difference in the plots being the 
values of the upper cut-offs due to the finite size of the 



As L — > cxD, the upper cutoff is oo because the maximum 
backbone mass is not constrained by the length of the 
shortest path. Substituting Eqs. (^), (|^), and (^ into 
Eq. (g), and equating powers of r (or powers of i) of the 
left and right hand sides of the resulting equation, we 
find 



ip ^ gt 



-{tb - 1). 



(6) 



Using Tb = rf/dg and the values gi = 2.04 |1^ and 
rfmin — 1-13 pi], we find ^|) = 1.72, in good agreement 
with our simulation result 



^ = 1.7 ±0.05. 



(7) 



V. AVERAGE CONDUCTANCE 

We can now calculate the average conductance. Since 
a is strongly correlated with (., and since a scales with r 
as r~^^ and i scales with r as r^'min, we have 

(cr) -^-'Vd,ni„. (8) 

Then using the fact that P{i\MB,r) ~ (.-'>' we have 

dl 



P{a\MB,r) - P{i = a-'^rnin/f^\MB, r) 

where 

z = (V-l)(d,„i„//i)-l--0.17. 
Now P{£\Mb, r) is nonzero only for 



dcr 



(9) 

(10) 

(11) 



where a and b are constants. Hence using (a) ^ i i^/'^min^ 
we find P{a\MB,r) is nonzero for 



(6Mb) 



('^min/''B)(-A/dmm) 



(&A/b) 



-/I/di 



(12) 



Using these bounds to normalize the distribution, we find 



P{a\MB,r) 
Then 

{a{MB,r))^ 



(ar)-A(2+i) _ ^fjM^Y-f./dsKz+i) ■ 



(13) 



Z+l 



aP{a\MB,r)da 



1- 



-(ar)- 



(bMB 



-A(z+2) 



1- 



(bMB 



-m(z+1) 



(14) 



Thus as Mb goes to infinity, {a{MB, r)) decreases asymp- 
totically to a constant as 



HMB,r)) 



z + l 



(ar)- 



{bAlB) 



l/dE 



-A(z+1)- 



(15) 



By considering the asymptotic dependence of {a{MB,r)) 
on Mb, we can reasonably fit the simulation results by 
choosing the parameters a and b in Eq. ( [l4|) to be 0.9 and 
6, respectively. Using these values for a and b, in Fig. 0(a) 
we plot (a) from Eq. ( |l4| ) for multiple values of r and find 
that agreement with the simulation results improves with 
increasing r. For large r, the curves for the simulation re- 
sults and the curves for the theoretical results are coinci- 
dent at large Mb- The poor results for small r are due to 



corrections-to-scaling not included in our derivation (e.g., 
for small r, there are significant corrections-to-scaling for 
the relations a - r^P^ and Mb ^ r'^^ ^^). 

Equation (|lj) can be re-cast in terms of the scaled 
variable x = MbJt^^ as 



z + l 



(16) 



where 



/(^) 



1 - (Jf-) 



-(p/dB)iz+2} 



1 - (^-) 



-{P./dB){z + l)' 



(17) 



In Fig. 0(b), we plot the average conductance scaled in 
accordance with Eqs. (^ and (p^). The expected col- 
lapse improves with increasing r for the same reason as 
noted above. 

Above the percolation threshold, for backbones of size 
larger than the correlation length, the strong correlation 
between the conductance and the shortest path breaks 
down and we expect the conductance to increase with 
the mass of the backbone, as is the case in non-random 
systems. This is seen in Fig. @(c), where we plot con- 
ductance versus backbone mass for the bond occupation 
probabilities p — 0.56 and p — 0.60, which are above the 
percolation threshold and, for comparison, conductance 
at the percolation threshold, p = 0.50 [^. Figure ||(c) 
shows that for p = 0.60, all backbone masses sampled 
are of size greater than the correlation length and the 
conductance increases monotonically. For p = 0.56, the 
smaller backbone masses are of size less than the corre- 
lation length and Fig. 0(c) shows that the conductance 
initially decreases; for larger backbone masses, however, 
the sizes of the backbones are greater than the corre- 
lation length and Fig. 0(c) shows that the conductance 
then increases. 



VI. DISCUSSION 



The derivation of Eq. ( |14[ ) and its agreement with the 
results of our simulations confirm our understanding of 
why the average conductance decreases with increasing 
backbone mass: the smaller contributions to the average 
conductance from the longer minimal paths possible in 
the clusters with larger backbone size cause the average 
conductance to be smaller. Our derivation was not spe- 
cific to two dimensions, and should also hold in higher 
dimensions. 
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Because, above the percolation threshold, the backbone 
mass is strongly correlated with the system size, we cre- 
ate backbones of a given mass by varying the system size. 
The simulations of Fig. 2(c) for the plots above the per- 
colation were obtained choosing L = 6, 10, 20, 40, 80, 
160, and 200. 
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(a) 

1=5 
(7-=0.20 



(b) 
M^=28 
1=5 
cr=0.22 
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M^=28 
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FIG. 1. Typical percolation clusters. The striped sites and the black sites are both part of the percolation cluster; only 
the black sites form the cluster backbone defined by sites A and B (black squares). Note in (a), with backbone size 6, the 
shortest path/, between A and B is 5 — the length of the shortest path must always be less than the backbone mass; (b) and 
(c) illustrate that in clusters with large backbones, the length of the shortest path between A and B can take on a broad range 
of values. 




FIG. 2. Average conductance versus backbone mass, (a) 
Simulation results at the percolation threshold for r = 8, 16, 
32, 64, and 128, where r is the distance between the two sites 
A and B\ the adjacent lines are the theoretical results. For 
large r, the curves for the simulation results and the corre- 
sponding curves for the theoretical results coincide for large 
Mb- (b) Plots of scaled backbone conductance for r = 8, 
16, 32, 64, and 128, scaled in accordance with Eqs. (|l^) and 
(p7|). The solid line is a plot of Eq. (119) with parameters a 
and h chosen as 0.9 and 6 respectively to best fit the values 
for r — 128 (right-pointing triangles). The collapse to this 
line for the lower values of r improves with increasing r. (c) 
Average conductance for r = 4 for p — 0.50, 0.56, and 0.60. 
For p = 0.56, the conductance as a function of Mb is not 
monotonia but rather has a minimum indicated by the arrow. 
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FIG. 3. (a)-(c) Three generations of the Sierpinski gasket. 
At each successive generation, there are two types of pairs of 
points: pairs which correspond to pairs in the previous gen- 
eration and pairs which do not. For example, in the second 
generation the pairs AB and AC are present in the previous 
generation but the pair DE has no corresponding pair in the 
previous generation. Similarly, the pairs AB, AC, and DE in 
the third generation correspond to pairs in the second gener- 
ation, but the pair EG does not. Because all points are mul- 
tiply connected, the mass of the backbone between any two 
points in each generation is equal to the mass of the entire 
gasket, (d) Average conductance between all pairs of points 
separated by distance r on the Sierpinski gasket versus the 
gasket mass. The points correspond to successive generations 
of the Sierpinski gasket. 
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FIG. 4. Distribution of shortest paths between A and B. 
(a) P{£\MB,r), the probability that the length of the short- 
est path between two points separated by distance r is £ for 
a given backbone mass, Mb- All plots are for r = 1, for 
values of Mb of 6 (squares), 96 (diamonds), 1536 (circles), 
and 24576 (triangles). The plots for the various Mb differ 
by the points at which they cut off. (b) P{£\Mb, r) for r = 1 
(squares), 4 (diamonds), and 16 (circles) for a single backbone 
size of 24576. (c) When scaled by r'*min, the plots collapse. 
The dashed line is constructed to have a slope of —1.7; see 
Eq. (a. 
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